CAFPE-123/09 
UG-FT-253/09 



Theory of the Hadronic Light-by-Light Contribution to 

Muon g-2^ 



Joaquim Prades 

CAFPE and Departamento de Fisica Teorica y del Cosmos, Universidad de 
Granada, Campus de Fuente Nueva, E-18002 Granada, Spain. 



Abstract 

I report on the theory, recent calculations and present status of the hadronic 
light-by-light contribution to muon g — 2. In particular, I report on work 
done together with Eduardo de Rafael and Arkady Vainshtein where we 
get a^^^^ = (10.5 ± 2.6) x 10"^'' as our present result for this quantity. 
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1 Introduction 



There are six possible momenta configurations contributing to the hadronic hght- 
by-hght to muon g-2, one of them is depicted in Fig. 1 and described by the vertex 
function 
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where g ^ is the momentum of the photon that couples to the external magnetic 
source, q = P2 — Pi = —ki — k2 — k^ and m is the muon mass. 




Figure 1: Hadronic light-by-light scattering contribution. 



The dominant contribution to the hadronic four-point function 

U'""'''iq,k,,k,,k2) = 

J dV J d^y J d^z e^i-'^^-+kry+k,-z) (^o\T[V^{0)V''{x)VP{y)V''{z)] |0) 

(2) 

comes from the three light-quark (g = u, d, s) components in the electromagnetic 
current V^^(x) = qQj'^q (x) where Q = diag(2, — 1, — 1)/3 denotes the light- 
quark electric charge matrix. For g — 2 we are interested in the limit g — where 
current conservation implies 
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Therefore, the muon anomaly can then be extracted as 

xtr{(^+m)[7/„7A](^+m)7j,(^+ h - "mY^lpil^- h-'ni)~^la} ■ 

(4) 

Here I report on the results of [1] and [2]. Previous work on the hadronic light- 
by-light contribution to muon g — 2 can be found in [3-12] and recent reviews are 
in [13-16]. 

The hadronic four-point function n^'^^'^(g, /ci, fca, /C2) is an extremely difficult 
object involving many scales and no full first principle calculation of it has been 
reported yet -even in the simpler large numbers of colors Nc of QCD hmit. No- 
tice that we need that hadronic four-point function with momenta /ci, k-2, and kz 
varying from to cxd and g — 0. Unfortunately, unlike the hadronic vacuum 
polarization, there is neither a direct connection of oF^^^ to a measurable quan- 
tity. Two lattice groups have started exploratory calculations [17, 18] but the 
final uncertainty that these calculations can reach is not clear yet. 

Attending to a combined large number of colors of QCD Nc and chiral per- 
turbation theory (CHPT) counting, one can distinguish four types of contribu- 
tions [19]. Notice that we use the CHPT counting only for organization of the 
contributions and refers to the lowest order term contributing in each case. In 
fact, Ref. [1] shows that there are chiral enhancement factors that demand more 
than Nambu-Goldstone bosons in the CHPT expansion in the light-by-light con- 
tribution to the muon anomaly. See more comments on this afterwards. 

The four different types of contributions mentioned above are: 

• Nambu-Goldstone boson exchanges contribution are 0{N^ and start con- 
tributing at 0{p^) in CHPT. 

• One-meson irreducible vertex contribution and non-Goldstone boson ex- 
changes contribute also at 0{N^ but start contributing at 0{p^) in CHPT. 

• One-loop of Goldstone bosons contribution are 0{\/N^ and start at 0{p'^) 
in CHPT. 

• One-loop of non-Goldstone boson contributions which are 0{1/Nc) but start 
contributing at in CHPT. 

Based on the counting above there are two full calculations [3,4,6] and [5,7]. There 
is also a detailed study of the 7r° exchange contribution [8] putting emphasis in 
obtaining analytical expressions for this part. Recently, two new calculations 
of the pion exchange using also the organization above have been made. In 
Ref. [10], the pion pole term exchange is evaluated within an effective chiral 
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model, NxQM. These authors also study the box diagram one-meson irreducible 
vertex contribution. The results are numerically very similar to the ones found in 
the literature as can be seen in Table 1. In Ref. [11], the author uses a large Nc 
model 7r°7*7* form factor with the pion also off-shell. This has to be considered 
as a first step and more work has to be done in order to have the full hght-by-light 
within this approach. In particular, it would be very interesting to calculate the 
contribution of one-meson irreducible vertex contribution within this model. 

There is also model independent short-distance QCD information on the rele- 
vant form factor. Using operator product expansion (OPE) in QCD, the authors 
of [12] pointed out a short-distance constraint of the reduced full four-point Green 
function (form factor) 

(0|T [V''{h)VP{h)V''{-{h + k2 + q))] |7(?)) (5) 

when q' — > and in the special momenta configuration —ksf ~ — A;| >> —{ki + 
k^Y Euchdean and large. In that kinematical region, 

T [V'ik^W^ik^)] ^ i e^P^^K \qQ^ lal, q] (6) 



k 

with A; = (/ci — /c3)/2 ~ /ci ~ —k^, . See also [20]. This short distance constraint 
was not explicitly imposed in calculations previous to [12]. 



2 Leading in 1/Nc Results 

Using effective field theory techniques, the authors of [9] shown that the leading 
large Nc contribution to a^^^^ contains an enhanced log {Mp/m.,j) term at low 
energy. Where the rho mass Mp acts as an ultraviolet scale and the pion mass m^r 
provides the infrared scale. The leading logarithm term is generated by Nambu- 
Goldstone boson exchange contributions and is fixed by the Wess-Zumino-Witten 
(WZW) vertex 77^77. 

™ = ©'^.S[ln^^..(ln^^).H ,r) 

In the chiral limit, where quark masses are neglected, and at large Nc-, the 
coefficient of this double logarithm is model independent and has been calculated 
and shown to be positive in [9] . All the calculations we discuss here agree with 
these leading behaviour and its coefficient including the sign. A global sign 
mistake in the tt^ exchange in the results presented in [3-5] was found by [8,9] 
and confirmed by [6,7] and by others [21,22]. The subleading ultraviolet scale 
/^-dependent terms [9], namely, log(/i/m7r) and a non-logarithmic term fi;(/x), are 
model dependent and calculations of them are implicit in the results presented 
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Figure 2: A generic meson exchange contribution to the hadronic hght-by-hght 
part of the muon g — 2. 



in [3-5,7,12]. In particular, ^(/i) contains the large Nc contributions from one- 
meson irreducible vertex and non-Nambu-Goldstone boson exchanges. In the 
next section we review the recent model calculations of the full leading in the 
1/Nc expansion contributions. 

2.1 Model Calculations 

The pseudo-scalar exchange is the dominant numerical contribution and was sat- 
urated in [3-8,10,11] by Nambu-Goldstone boson exchange. This contribution 
is depicted in Fig. 2 with M = vr^,?], r/'. The relevant four-point function was 
obtained in terms of the off-shell (A;i)7*(A;3) form factor J^{k\,k1) and the 
off-shell 7r'^7*(fc2)7('? = 0) form factor JF(A;2,0) modulating each one of the two 
WZW 7r°77 vertex. 

In all cases discussed here, several short-distance QCD constraints were im- 
posed on these form-factors. In particular, they all have the correct QCD short- 
distance behaviour 

HQ\Q')^^, and ^(g^o)^^ (8) 

when Q"^ is Euclidean and large and are in agreement with 7r°7*7 low-energy data 
^. They differ slightly in shape due to the different model assumptions (VMD, 
ENJL, Large Nc-, NxQM) but they produce small numerical differences always 
compatible within quoted uncertainty ~ (1.3 — 1.6) x 10~^° -see Table 1. 

^See however the new measurement of the 77* — > ttq transition form factor by BaBar [23] 
at momentum transfer energies between 4 GeV^ and 40 GeV^ 
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References 


10 


10 X a 




7r° only 


7r°, T) and r]' 


[3,4,6] 


5.7 


8.3 ± 0.6 


[5,7] 


5.6 


8.5 ± 1.3 


[8] with /i2 = 


5.8 


8.3 ± 1.2 


[8] with /i2 = -10 GeV^ 


6.3 




[10] 


6.3 - 6.7 




[11] 


7.2 


9.9 ± 1.6 


[12] 


7.65 


11.4±1.0 



Table 1: Results for the tt", rj and rj' exchange contributions. 



A [GeV] 


0.7 


1.0 


2.0 


4.0 




2.2 


2.0 


1.9 


2.0 



Table 2: Sum of the short- and long-distance quark loop contributions [5] as a 
function of the matching scale A. 



Within the models used in [3-8, 10, 11], to get the full contribution at leading 
in 1/Nc one needs to add the one-meson irreducible vertex contribution and the 
non-Goldstone boson exchanges. In particular, below some hadronic scale A, 
the one-mcson irreducible vertex contribution was identified in [5, 7] with the 
ENJL quark box contribution with four dressed photon legs. While to mimic the 
contribution of short- distance QCD quarks above A, a loop of bare massive heavy 
quark with mass A and QCD vertices was used. The results are in Table 2. There, 
one can see a very nice stabihty region when A is in the interval [0.7, 4.0] GeV. 
Similar results for a constituent quark-box contribution below A were obtained 
in [3,4], though these authors didn't discuss any short-distance-long-distance 
matching. 

In [5,7], non-Goldstone boson exchanges were saturated by the hadrons ap- 
pearing in the model, i.e. the lowest scalar and pseudo-vector hadrons. There, 
both states were used in nonet-symmetry -this symmetry is exact in the large 
limit of QCD. 

Within the ENJL model, the one-meson irreducible vertex contribution is 
related trough Ward identities to the scalar exchange which we discuss below 
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References 


IQio X a™ 


[3,4,6] 


0.17 ± 0.10 


[5,7] 


0.25 ± 0.10 



Table 3: Results for the axial- vector exchange contributions from [3,4,6] and [5,7]. 




Point — Like 



Figure 3: Goldstone boson exchange in the model in [12] contributing to the 
hadronic light-by-light. 



and both have to be included within this model [5,7]. The result of the scalar 
exchange obtained in [5] is 

a"^^^(Scalar) = -(0.7 ± 0.2) x 10"^° . (9) 

The scalar exchange was not included in [3,4,6,8]. The result of the axial- vector 
exchanges in [3,4,6] and [5,7] can be found in Table 3. 

Melnikov and Vainshtein used a model that saturates the hadronic four-point 
function in (2) at leading order in the 1/Nc expansion by the exchange of the 
Nambu-Goldstone TT^^rj^rj' and the lowest axial-vector fi states. In that model, 
the new OPE constraint of the reduced four-point function found in [12] men- 
tioned above, forces the T^^^*[q)^{pz = 0) vertex to be point-like rather than 
including a jF(g^,0) form factor. There are also OPE constraints for other mo- 
menta regions [24] which are not satisfied by the model in [12] though the authors 
argued that this mismatch makes only a small numerical difference of the order 
of 0.05 X 10~^°. In fact, within the large Nc framework, it has been shown [25] 
that in general for other than two-point functions, to satisfy fully the QCD short- 
distance properties requires the inclusion of an infinite number of narrow states. 
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3 Next-to- leading in l/Nc Results 



For the next-to-leading in 1/Nc contributions to the a^^^^ there is no model inde- 
pendent result at present and is possibly the most difficult component. Charged 
pion and kaon loops saturated this contribution in [3-7]. To dress the photon 
interacting with pions, a particular Hidden Gauge Symmetry (HGS) model was 
used in [3,4,6] while a full VMD was used in [5,7]. The results obtained in these 
two models are -(0.45 ±0.85) x lO^^" in [3] and -(1.9 ±0.5) x 10-^° in [5] while 
using a point-like vertex one gets —4.6 x 10^^°. 

Both models (HGS and VMD) satisfy the known constraints though start 
differing at 0{p^) in CHPT. Some studies of the cut-off dependence of the pion 
loop using the full VMD model was done in [5] and showed that their final number 
comes from fairly low energies where the model dependence should be smaller. 

The authors of [12] analyzed the model used in [3,4] and showed that there 
is a large cancellation between the first three terms of an expansion in powers 
of {mTr/MpY and with large higher order corrections when expanded in CHPT 
orders but the same applies to the 7r° exchange as can be seen from Table 6 in 
the ffist reference in [2] by comparing the WZW column with the others. The 
authors of [12] took (0 ± 1) x 10^^° as a guess estimate of the total NLO in 
1/Nc contribution. This seems too simply and certainly with underestimated 
uncertainty. 

4 Comparing Different Calculations 

The comparison of individual contributions in [3-8, 10-12] has to be done with 
care because they come from different model assumptions to construct the full 
relevant four-point function. In fact, the authors of [10] have shown that their 
constituent quark-box provides the correct asymptotics and in particular the new 
OPE found in [12]. It has more sense to compare results for a^^^^ either at leading 
order or at ncxt-to-leading order in the l/N^ expansion. 

The results for the final hadronic light-by-light contribution to a^^^^ quoted 
in [3-7,12] are in Table 4. The apparent agreement between [3,4,6] and [5,7] 
hides non-negligible differences which numerically almost compensate between 
the quark-loop and charged pion and [12] are in Table 4. Notice also that [3,4,6] 
didn't include the scalar exchange. 

Comparing the results of [5,7] and [12], as discussed above, we have found 
several differences of order 1.5 x 10~^° which are not related to the new short- 
distance constraint used in [12]. The different axial- vector mass mixing accounts 
for —1.5 X 10~^°, the absence of the scalar exchange in [12] accounts for —0.7 x 
10-1° and the use of a vanishing NLO in l/N^ contribution in [12] accounts for 
—1.9 X 10"^°. These model dependent differences add up to —4.1 x 10"^° out of 
the final —5.3 x 10"^° difference between the results in [5,7] and the ones in [12] 
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Full Hadronic Light-by-Light 


10^° X 


[3,4,6] 


8.9± 1.7 


[5,7] 


8.9 ± 3.2 


[12] 


13.6 ± 2.5 



Table 4: Results for the full hadronic light- by-light contribution to a^^^^. 

-see Table 4. Clearly, the new OPE constraint used in [12] accounts only for a 
small part of the large numerical final difference. 

5 Conclusions and Prospects 

To give a result at present for the hadronic light-by-light contribution to the 
muon anomalous magnetic moment, the authors of [1] concluded, from the above 
considerations, that it is fair to proceed as follows: 
Contribution to a^^^^ from 7r°, r) and r)' exchanges 

Because of the effect of the OPE constraint discussed above, we suggested [1] to 
take as central value the result of Ref. [12] with, however, the largest error quoted 
in Refs. [5,7]: 

a"L'^^(7r ,ri,ri') = (11.4 ± 1.3) x IQ-^" . (10) 

Recall that this central value is quite close to the one in the ENJL model which 
includes the short-distance quark-loop contribution. 
Contribution to a^^^^ from pseudo-vector exchanges 

The analysis made in Ref. [12] suggests that the errors in the first and second 
entries of Table 3 are likely to be underestimates. Raising their ±0.10 errors to 
±1 puts the three numbers in agreement within one sigma. We suggested [1] then 
as the best estimate for this contribution at present 

a"^''^(pseudo - vectors) = (1.5 ± 1) x 10~^° . (11) 

Contribution to a^^^^ from scalar exchanges 

The ENJL-model should give a good estimate for these contributions. We kept 
[1], therefore, the result of Ref. [5,7] with, however, a larger error which covers 
the effect of other unaccounted meson exchanges, 

a"^'^^(scalars) = -(0.7 ± 0.7) x 10"^° . (12) 

Contribution to a^^^^ from dressed charged pion and kaon loop 

Because of the instability of the results for the charged pion loop and unaccounted 
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loops of other mesons, we suggested [1] using the central value of the ENJL result 
but wit a larger error: 

a"^''^(7r-dressed loop) = -(1.9 ± 1.9) x 10^^° . (13) 

Prom these considerations, adding the errors in quadrature, as well as the 
small charm contribution 0.23 x 10~^°, we get 

^HLbL ^ 5 ^ 2.6) X 10"^° , (14) 

as our final estimate. 

The proposed new muon g — 2 experiments at Fermilab [28] with 1.6 x 10~^° 
accuracy goal and at J-PARC [29] with even higher accuracy goal between 1.2 x 
10"^'' and 0.6 x 10~^° call for a considerable improvement in the present cal- 
culations of a^^^^. The use of further theoretical and experimental constraints 
could result in reaching such accuracy soon enough. In particular, imposing as 
many as possible short-distance QCD constraints [3-8,11] has result in a better 
understanding of the numerically dominant 7C^ exchange. At present, none of 
the light-by-light hadronic parametrization satisfy fully all short distance QCD 
constraints. In particular, this requires the inclusion of infinite number of narrow 
states for other than two-point functions and two-point functions with soft in- 
sertions [25] . A numerical dominance of certain momenta configuration can help 
to minimize the effects of short distance QCD constraints not satisfied, as in the 
model in [12]. 

More experimental information on the decays 7r° — > 77*, 7r° — > 7*7* and 
7r° — > e+e" (with radiative corrections included [22,26,27]) in the low- and 
intermediate-energy regions (below a few GcVs) can also help to confirm some of 
the neutral pion exchange results. A better understanding of other smaller contri- 
butions but with comparable uncertainties needs both more theoretical work and 
experimental information. This refers in particular to pseudo-vector exchanges. 
Experimental data on radiative decays and two-photon production of these and 
other C-even resonances can be useful in that respect. 

New approaches to the pion dressed loop contribution, together with exper- 
imental information on the vertex 7r"'"7r~7*7* in the intermediate energy region 
(0.5 — 1.5) GeV would also be very welcome. Measurements of two- photon pro- 
cesses like e"^e~ e+e-Tr+Tr" can be useful to give information on that vertex 
and again could reduce the model dependence. The two-gamma physics program 
low energy facilities like the experiment KLOE-2 at DA$NE will be very useful 
and well suited in the processes mentioned above which information can help to 
decrease the present model dependence of a^^^^. 
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